A = str2sym('[1/x,1/(x+1);1/(x+2),1/(x+3)]');
B = str2sym('[x,1;x+2,0]');
C = B - A
D = A \ B
